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We study macroscopic quantum coherence (MQC) in small 
magnetic particles where the magnetization (in ferromagnets) 
or the Neel vector (in antiferromagnets) can tunnel between 
energy minima. We consider here the more general case of 
MQC in ferrimagnets by studying a model for a mesoscopic 
antiferromagnet with an uncompensated magnetic moment. 
Through semi-classical calculations we show that even a small 
moment has a drastic effect on MQC. In particular, there is 
a rapid crossover to a regime where the MQC tunnel split- 
ting is equal to that obtained for a ferromagnet, even though 
the system is still an antiferromagnet for all other aspects. 
We calculate this tunnel splitting via instanton methods and 
compare it with numerical evaluations. As an application we 
re-examine the experimental evidence for MQC in ferritin and 
show that even though the uncompensated moment of ferritin 
is small it greatly modifies the MQC behavior. The excess 
spin allows us to extract values for experimental parameters 
without making any assumption about the classical attempt 
frequency, in contrast to previous fits. Finally, we also discuss 
the implications of our results for MQC in molecular magnets. 
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I. INTRODUCTION 

Quantum phenomena, like tunneling, arc usually lim- 
ited to the atomic world. However, under some stringent 
conditions it is also possible for a macroscopic variable 
to exhibit quantum behavior. Over the last 15 years, 
a lot of experimental and theoretical work has been de- 
voted to the quest of such systems. The clearest evidence 
so far has been observed in Josephson junctions, jSthere 
quantum tunneling of the current phase can occuran 

More recently, another class of systems where macro- 
scopic quantum phenomena might b& .possible has been 
considered: small magnetic particlesB'u. In a grain too 
small to accommodate more than one magnetic domain, 
all the spins are locked together by the exchange inter- 
action, and only their global orientation can change. Be- 
cause of magnetic anisotropics and external fields, the 
grain energy depends on this orientation. Two type of 
phenomena can be envisaged. First, macroscopic quan- 
tum tunneling (MQT), where the magnetization tunnels 
through an energy barrier from one local minimum to 
the global one. Second, macroscopic quantum coherence 
(MQC), where the classical energy has several equiva- 
lent minima, and the magnetization oscillates back and 



forth between them. This oscillation should show up in 
magnetic noise and susceptibility spectraQ. 

Theoretical calculations based on simple models were 
first performed for ferromagnets (FM)ffiju. Surprisingly, 
they revealed that MQC was possible in systems con- 
taining as much as 10 5 moments. This generated a lot 
of interest. Later extensionsolid of these calculations to 
antiferromagnets (AFM), where the Neel vector is the 
tunneling entity, showed that MQC should show up at 
higher temperatures and higher frequencies than in FM 
grains of similar size. This makes AFM more interesting 
for experimental purposes. n| _ . 

The models used in early calculationsallj for AFM are 
based on the assumption that the two sublattices had 
an exactly equal number of sites, meaning that the grain 
had no net magnetic moment. However, this condition is 
not met experimentally in finite-size magnets, since be- 
cause of statistical fluctuations a small imbalance in thg 
two sublattices is expected. This was recognized earlyO 
and deemed as a practical advantage. As the uncom- 
pensated spin is locked together with the Neel vector, it 
should mirror the MQC oscillations; the resulting mag- 
netization oscillations are easier to detect than those p£ 
the Neel vector. Although addressed theoreticallyllirE£| 
the influence of this uncompensated moment on MQC 
for a physical system has not been studied yet, which 
motivated the present investigation. 

To get an idea about the magnitude of this uncompen- 
sated moment, let us consider several systems of prac- 
tical interest. Ferritin is a macromolecule of high bio- 
logical iiapoptaace, for which MQC evidence has been 
reportecSia'Ej 113. It has a magnetic core containing up 
to 4500 Fe 3+ atoms, whose spins are antiferromagnet- 
ically coupled. Yet a ferritin molecule at full loading 
has a magnetic moment corresponding to that of roughly 
50 iron atoms. Clearly, this is small on the scale of 
the particle. Other systems of interest are molecular 
magnetsEao. These molecules contain ligands binding a 
number of magnetic atoms (so far, up to 20), whose spins 
are strongly coupled. In addition, macroscopic samples 
of perfectly identical molecules can be prepared. This 
makes them interesting for MQC, because the tunnel- 
ing frequency scales exponentially with the number of 
atoms. In these molecules, competing ferro- and anti- 
ferromagnetic interactions usually result in ferrimagnetic 
structures with uncompensated moments that can be as 
large as, or even larger than, the sublattice magnetiza- 
tion. 

In this paper we will investigate the influence of the 
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uncompensated moment on MQC in ferrimagnets. In 
the two systems we have just discussed, the magnitude 
of this moment ranges from large (most molecular mag- 
nets) to seemingly negligible (ferritin). Nevertheless, we 
will show that even in ferritin uncompensated spins dras- 
tically alter the MQC results and dominate the tunneling 
dynamics. The paper is structured in the following way. 
In Section ||, we introduce our model for a ferrimagnet: 
an antiferromagnet coupled to an uncompensated mo- 
ment. We discuss its energy spectrum semi-classically by 
looking at its instanton solutions in the regime of inter- 
est for experiments. In Section HI, we apply our results 
to ferritin and molecular magnets. We conclude with a 
discussion in Section IV. 



(t) = ±2 arctan(exp(wAFMT)), 0{t) = it/2, (2) 



II. MODEL AND SEMI-CLASSICAL RESULTS 

In an antiferromagnetic grain much smaller than the 
typical domain wall width, the Neel vector is expected 
to be uniform (Stoner-Wolfarth model). Its low-energy 
dynamics is describedrbjVnp*i effective action analogous tei 
that of a rigid rotatorE9'BEil. This model was extendcdtil 
to include the effects of an uncompensated moment, 
which is treated as one (large) single spin coupled to the 
AFM core. In the semi-classical regime, aiid-|Using the 
coherent-state path-integral representationl3E2l, one gets 
the following (Euclidean) action in the limit whpce the 
excess spin is strongly coupled to the Neel vectoiO 



2 1 j> 2 sin 2 9) + K y sin 2 6 sin 2 . 



2 7 2 

+ K z cos 2 9^ +ihS I dr<P(l~ cos 9), (1) 



where 9 and <j> are the spherical coordinates of the Neel 
vector, V is the volume of the grain, x± is its perpendic- 
ular magnetic susceptibility, K z > K y > its magnetic 
anisotropiesa, 7 = 2 /ig /h, and S is the magnitude of 
the excess spin. We have therefore included an easy axis 
anisotropy along the s-axis, with the xy-plane taken as 
the easy plane. 

Notice that the term ih,S Jdr 4> in the action only de- 
pends on the initial and final states, and hence does not 
affect the classical equations of motion. Howeverliio, 
it associates a different phase factor exp(±2i7rS') to in- 
stantons and anti-instantons. This makes the semi- 
classical predictions conform with Kramers theorem: If S 
is half-integer, the contributions of instantons and anti- 
instantons exactly cancel each other, and the tunnel split- 
ting disappears. If S is integer, the contributions sum up 
coherently, and tunneling occurs. We will suppress this 
topological term in the following discussion of the classi- 
cal solutions but will reinstate it in the final expression 
for the tunnel splitting. 

In the absence of an excess spin, i.e. for S = 0, the 
instanton and anti-instanton solutions are given by 
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These solutions have an action 



SafmA = V^2K v x±/ Hi 



(3) 
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We-will be interested in the regime K z ^> K y . Loss et 
aZ.EII have investigated the effect of a small excess spin, 
small in the sense that HS/Safm <C K z /K y . Their 
results are based on approximations that can be self- 
consistently checked in this regime. They found that the 
preceding expressions for wafm and Safm remain correct 
provided one replaces x± by an effective susceptibility 
X f = X ± + (2/K z )(S f i B /V) 2 . 

se results were later confirmed by Duan and 
i3, using an alternative model for the uncompen- 
sated moment. Unlike Loss et at, Duan and Garg use 
a variational approach, which gives results for the entire 
parameter space. The drawback of the variational appo- 
rach is that it is not possible to assess the accuracy of its 
predictions, this being dependent on the ansatz quality. 

We will now discuss the instanton solutions of the 
action (|]J) for a large excess spin, more precisely for 
hS/S afm ^ 1- We will perform a number of approxima- 
tions, and will check their self-consistency at the end. 

In the absence of an excess spin, the characteristic time 
scale is w^J M . To simplify the equations of motion, we 
make a change of variables to a reduced time f = ijjafmt. 
To get finite-action solutions, (f> must be real, and 9 — n/2 
must be purely imaginary. Defining ii) = 9 — n/2, one 
gets for the action 

Se = ^± j d f f {-tf 2 + cosh 2 d{j> 2 + sin 2 4>) (5) 



^ sinh 2 ® - A- 
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where overdots now mean derivatives with respect to f . 
Notice that all the physics is determined by only two di- 
mensionless quantities, hS/S 'afm and K z /K y . The equa- 
tions of motion are 



i sin 2(f) - 



ns 



cosh ■& Safm 



i?-2tf<£tanhtf 



(7) 



K z 



4 = ^ sinh 2tf [~4> 2 + — - sin 2 , 



Pi Q 

2- ^coshtf, (8) 

<~>AFM 



with the boundary conditions for an instanton solution 
being <p(— 00) = 0, 4>(+oo) = ir, #(±00) = 0. 

Recalling that we are interested in the regime 
hS/S afm > 1, K z /K y » 1 we see that in Eqs. (0) and 
(||) the second-order time-derivatives have much smaller 
prefactors than the first-order ones. Therefore, we ex- 
pect the second-order derivatives to be negligible, except 
in possible boundary layers, small regions usually located 
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close to the boundaries (though not always, see the dis- 
cussion in Ref. It turns out, however, that there 
are no boundary layers for our choice of boundary con- 
ditions. If we neglect the i? term in Eq. (||), and assume 
that cj> 2 + sin 2 cf> -C K z /K v , we get 

= -arcsinh ( 2^-^-^) . (9) 
V K z Aafm / 

We will see later that {K y / K z ){TiS/ Safm)4> < 1. We can 
therefore expand the hyperbolic arcsinus to first order. 
Substituting then ■& into Eq. (0), neglecting c/>sinh$ with 
respect to TiS/Safm, and expanding cosh?9 to leading 
order, we obtain 

( 1+4 t(^) 2 )* = ^'°* (10 » 

Thus, one gets exactly the same equation for 4>{ T ) as m 
the absence of the excess moment, except that now the 
instanton frequency is given by 

-( i+ £(&)T 

We see that the uncompensated moment slows down the 
transition between the two energy minima (by a fac- 
tor given precisely by u>). It is now straightforward to 
check-ilie self-consistency of the approximations that we 
usecM 

Moreover, it is easy to check that, within the same level 
of approximation, w also gives the frequency of the small 
oscillations about each minimum and thus describes the 
lowest energy level spacing (in the absence of the tunnel- 
ing splitting) in the wells. 

If we reinstate the full units, the instanton frequency 
in the ferrimagnetic regime (S ^ 0) becomes 



WPerri = -j^-sj K y K z , (12) 

with 

For the calculation of the action, sinceEl $ is at most of 
order */ K y /K z , we keep only the lowest-order terms in 
the Taylor expansions of sinhz? and coshi?. This gives 

Let us discuss these results in several limits. First, for an 
excess spin becoming large, i.e. for S — > 00, one sees that 
A — > 1. Hence, WFcm and S'Fcrrj^gpd to the expressions 
expected for pure ferromagnetM'W (with K z ^> K y ), 



2Y / 

wfm = j^y/K y K z , S F M = 2hSyK y /K z . (15) 

However, the role of the total spin in the ferromagnetic 
case is now played by the excess spin S\ To recover the 
FM regime starting from an AFM might seem at first 
surprising, but should have been expected. As the im- 
balance between the two sublattices grows, the dynamics 
becomes that of a large spin coupled to a small AFM, 
i.e. a ferromagnetic dynamics. This continuous transition 
between the seemingly unrelated FM and AFM regimes 
provides a useful check of the model. 

In the opposite limit of a small excess spin, i.e. for 
S — > 0, we would not expect our results to be accu- 
rate, since they were derived under the assumption that 
HS/S 'afm 3> 1. However, it is easy to check that our 
expression for WFcrri is in fact equal to that derived in 
Ref. [ll], where x± m Eq. ® 1S replaced by an effective 
susceptibility X f = X± + {2/K z ){Sfi B /V) 2 . 

Therefore, our results smoothly connect to the small 
excess spin regime. Our analytical expressions are valid 
over the whole range of S (but for K z ^> K v ), i.e. all the 
way from the AFM to the FM regime. An alternative 
way to interpret them is to notice that in our expression 
for ^Ferri, the term A (Safm/^-S) 2 lies between and 
K y /K Zl and hence is a small correction. The remaining 
term is exactly that expected for a pure ferromagnet, but 
with an effective spin S/ A; in a sense, all MQC splittings 
are given by the FM results. 

The various regimes can be characterized by the 
value of A. For a small excess spin, such that A ss 
2 ^K y /K z (hS/S afm) ~ 0, we are in an AFM regime 
where (to first order) only \± an d K y enter the action 
and the instanton frequency. For a large excess spin, such 
that A ~ 1, we are in a FM regime, where (again to first 
order) only K y and K z are relevant. The crossover oc- 
curs for TiS/S afm ~ \jK z jK y . We can get an estimate 
of the critical value of the uncompensated moment by 
noting that, if MQC is to be observable, Safm/Ti should 
at most be of the order of 10. For a ratio K z /K y of about 
10, the crossover would occur for S ~ 30, i.e. an uncom- 
pensated moment of 60/i_b- This excess moment is quite 
small and a large enough system might be indistinguish- 
able from a pure AFM for all its properties except for its 
MQC behavior. 

To assess the accuracy of our analytical predictions, 
and in particular to estimate how large the ratio K z /K y 
needs to be, we can compare them with direct numerical 
calculations of the instanton solution of Eqs. (0) and (||). 
Let us first fix the ratio K z /K y to a value of 30, and vary. 
HS/S afm- We plot in Fig. [l] the instanton frequencyEll 
and action. As expected in that regime, the analytical 
and numerical results compare extremely well. Trying 
now with a much smaller ratio K z /K y = 5, where our 
approximations are not so well justified, we get the re- 
sults plotted in Fig. ^. Our analytical results for the 
instanton frequency are accurate for small and large val- 
ues of HS/S afm', deviations are apparent (and expected) 
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in the intermediate regime. For the action, good results 
are obtained for small HS/Safm] for large values the dis- 
crepancies are due to our first-order expansion of sinh6* 
and cosh# in the action. We have verified that inclusion 
of higher-order terms gives excellent agreement between 
analytical and numerical results. 

Alternatively, we can fix HS/Safm and vary K z /K y . 
The results for the choice HS/Safm = 10 are shown in 
Fig. [5| Even though our analytical expressions were de- 
rived under the assumption that K z K y , we see that 
they give accurate predictions down to K z /K y ks 2. For 
still lower values the true action diverges, whereas our 
expressions remain finite. This divergence reflects the ap- 
pearance of a new symmetry: For K y — K z , S'x is 9, con- 
served quantity, and consequently tunneling is blocked if 
5^0. The tunnel splitting must therefore vanish. 

With our solution for the instanton, and after evalua- 
tion of the fluctuation determinant (again for K z ^> K y ) 1 
we get the following prediction for the tunnel splitting 



8hu 



Fcrri 



"JFerri 
2nh 



cos{irS)\e- SF °" i/n . 



(16) 



The term |cos(7r5)| is duc0 to the topological term, 
which we have reinstated. 

Using the criterion that the MQC resonance should 
show up at temperatures low enough that the quan- 
tum transition rates between the degenerate energy min- 
ima are higher than the classical rates, one gets for the 
crossover temperature 



k B T* = K y VK/S FeIli . 



(17) 



Up to relative corrections of order K y /K z , Sperri is equal 
to TKyF/wFerri, and therefore k B T* w Since 
the two lowest levels are separated from the others by an 
energy Tk^pcrrii the criterion means that the transition to 
the quantum regime occurs for temperatures such that 
only the lowest two levels are populated. 

To conclude this section, let us come back to the issue 
of the best system to detect MQC, ferromagnets or anti- 
ferromagnets. Early calculations of the tunnel splitting in 
(pure) FM and AFM gave A FM ~ uj exp(-Nsy/ K y / K z ) 
and A afm ~ ^>o exp(— Ns^fKy/ J), where N is the num- 
ber of sites, s is the magnitude of the spins, J is the ex- 
change coupling, and loq a characteristic frequency. Since 
J ^> K z , MQC should be easier to observe in AFM than 
in FM. Our results for an AFM with an uncompensated 
moment give Apem ~ exp{—N e Sy/K y /K z ), where N e 
is the number of excess spins. Even though this is similar 
in structure to the FM result, AFM nevertheless remain 
more favorable to observe MQC since typically N e <C N. 



III. APPLICATIONS 



A. Ferritin 

Ferritin is a macromolecule of high biological im- 
portance, naturally occurring in many living organisms 
where it is responsible for iron transport. It has a chest- 
nut structure, with an outside shell about 2nm thick 
enclosing a core of iron oxide. The roughly spherical 
core has a diameter of 7.5 nm, and can contain a variable 
number .-of irons, up to a maximum of 4500.^-iThie oxide 
appearsca structurally similar to ferrihydritccjEd . 

Ferritin moJfiGules have a small magnetic moment. 
Measurements&E3 on fully loaded grains of horse-spleen 
ferritin give values ranging between 217/is and 316/i B . 
This range corresponds to roughly 45-60 parallel iron 
spins, and, if compared to the moment expected if all 
the spins were parallel (about 4500 x 5/is), provides ev- 
idence for an antiferromagnetic core. .. .. . 

MQC evidence has been reportedlMfeJiaia in di- 
luted samples of horse-spleen ferritin. Since the reso- 
nance frequency depends exponentially on the amount 
of AFM material, it is crucial to be able to select par- 
ticles with a welLdefiaed loading. For that reason, the 
first experimentsQ'llilO used fully loaded grains, which 
could be reliably sifted out. In both the magnetic noise 
and susceptibility spectra, a peak at v rca — 940 kHz ap- 
peared below T* = 200 mK, and was fully developed 
at about 30 mK. From its field- and concentration- 
dependence, and from other cross-checks, this resjmance 
was attributed to MQCQ. Later experimentsEjiij ex- 
ploited recent advances in synthesis of ferritin to arti- 
ficially engineer grains with a given loading. The res- 
onance frequency was found to decrease exponentially 
with the loading, as expected from the MQC predictions 
for a (pure) AFM. 

The prftY. iqus . . quan titative analysis of the experimen- 
tal resultsBliJliil'Olia used the following estimate for the 
MQC resonance frequency, a; res = wo exp(— Safm/^), 
with ojq fixed to a value of 10 10 s _1 (and not measured or 
calculated from the inferred x± and K y ). The crossover 
temperature was estimated as the one where the rate of 
thermal fluctuations between the two minima is equal to 
the MQC rate, giving! T* = fi B ^/K v /2 X ±/k B . From 
the experimental results for u> Ies and T* , the following 
values were derived 



X-i 



5.2 ■ 10" 



emu 
Gem 3 



0.95 -10 3 -^Ir. (18) 
cm 3 



An assumed value for uio had to be used for lack of ex- 
perimental data at that time. Notice that uo/^tt gives 
the frequency of the small oscillations around the min- 
ima. In Kramers approach, it should therefore approxi- 
mately equal tbs-elassical attempt frequency fo. Recent 
measurementsSaS of /o in fully loaded grains of horse- 
spleen ferritin gave much larger values, ranging between 
10 11 and 10 12 s _1 . The spread in these results is indica- 
tive of the difficulty of these measurements, but they both 
point to values of ujq much larger than 10 10 s~ 1 . 
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This would hardly matter if MQC measurements had 
been performed on fully loaded grains only. Slightly dif- 
ferent values for the experimental parameters \± an d K y 
would result from the fit. However, uq plays an impor- 
tant role in the scaling of the resonance frequency with 
the grain loading, and within the pure AFM approach, 
where Lu rcs — luq exp(— Safm/K), it is readily seen that 
a value of 10 10 s _1 is necessary to reproduce the experi- 
mental data. Since this value is fat-too small compared 
with the independently measurecOLj one, one is there- 
fore faced with a consistency problem, which we shall 
address now. 

Being so small on the scale of the molecule, the uncom- 
pensated moment was neglected in the previous quanti- 
tative analysis. That it should play an important role 
can be easily appreciated. If one takes at face value 
the parameters (1^), one gets for the action a value 
Safm/Ti ~ 7.5. With a measured excess spin So ~ 100- 
150, the ratio US/Safm would be in the range 13-20. 
Thus, it is not possible to neglect the uncompensated 
moment. 

However, the following problem now arises. Experi- 
mentally, two quantities were measured in the MQC ex- 
periments: the resonance frequency ^ res and the crossover 
temperature T*. Without excess spin, these quantities 
depend on two parameters only (K y and Xx), which 
are therefore unambiguously fixed. With an excess spin, 
three parameters are relevant (K y , K z , and x±)i and 
the MQC experimental results alone are not sufficient to 
determine them. 

We can circumvent this problem by using an additional 
parameter, namely .the, Neel temperature Tjy. We use 
the simple estimateerE-a x± = where N 

is the number of spins (about 4500 for a fully loaded 
grain), and V is the core volume (which we estimate as 
that of a sphere of 7.5 nm in diameter). From the value 
Tn = 240 K reported for ferritir£3, one gets XJ- = 5 ' 
10~ 5 emu/Gem 3 . 

If we suppose that the excess spin is So = 100, Eqs. 
@ and (0) give 



q erg 
K v = 1.5- 10 3 — %. 

enr 



K z = 9.8 • 10 



5 er § 



(19) 



The value for K y is a bit larger than the one obtained 
neglecting the uncompensated moment; K z was not de- 
termined in the original fit. These results depend on the 
values taken for So and x±- For example, repeating the 
fits with S = 150 we find that K z = 2.2 • 10 6 erg/cm 3 , 
whereas with x± — 10~ 5 emu/Gem 3 , K z reduces to 
4.7 • 10 5 erg/cm 3 , while K v is left unchanged. 

From the parameters (|l9[), one gets for a fully loaded 
grain WFerri = 10 1:L s _1 , SVcrri ~ 12?i, and A w 0.65. We 
are therefore not quite in the FM regime, which explains 
why the choice of x± still has an influence. The value. p£ 
WFerri/27r is a bit lower than the experimental resultscSH 
for the classical attempt frequency fo (between 10 11 and 
10 12 s _1 ), but not unreasonably so considering the sim- 
plicity of the model and the experimental uncertainties. 



To get predictions for the volume dependence of the 
MQC resonance frequency, we need first to know how the 
excess spin depends on the grain loading. To the best of 
our knowledge, no experimental study exists. We will 
postulate a simple scaling law, S — So(V/Vo) u , where 
S and V are the excess spin and volume of a grain, So 
and Vq corresponding to a fully loaded grain. For the 
exponent v, several models have been proposed. Suppose 
the uncompensated moment is solely due to a difference 
in the number of-A and B sites at the surface. Barbara 
and Chudnovskycl have argued that, since the number of 
surface spins scales like V 2 ^ 3 , the excess spin (being the 
difference of two such random variables) should therefore 
scale like V 1 ' 3 . Neel has proposed several other modekfa, 
giving a scaling of S with either V 1 / 3 , V 1 / 2 , or V 2 / 3 . The 
theoretical picture is therefore far from settled and in the 
absence of experimental data we will consider all three 
scaling laws. 

The predictions for the MQC frequency as a function of 
the number of spins are shown in Fig. [|, together with the 
experimental measurements. A scaling of S with y 1 / 3 
clearly gives the best results, whereas a scaling with I/ 1 / 2 
predicts frequencies that are too large at low loadings. 
As this gets even worse with a F 2 / 3 -scaling, we have not 
displayed these results in Fig. || to avoid cluttering. 

Notice that the scaling is more complicated than in 
the simple AFM picture used in the previous analy- 
sis of the experiments. In the latter 
ujo exp(— Safm/K), with Safm being linear in the total 
number of spins. Hence, the resonance frequency de- 
pends exponentially on the loading. This is no more the 
case in the ferrimagnetic case. In the FM (A « 1) and 
AFM (A « 0) limits, the frequency does indeed depend 
exponentially on the loading, though in a different way 
in the two regimes. In the intermediate regime, there is 
a smooth crossover between the two scalings. 

The value taken for the uncompensated spin So has 
nearly no effect on the scaling predictions. (We did not 
show in Fig. [| the results obtained with So = 150, the 
curves being indistinguishable.) On the other hand, the 
value taken for x± affects the scaling markedly. We 
show in Fig. the predictions obtained with x± — 
10 -5 emu/Gem 3 , which suggests an alternative ap- 
proach. We recall that we have estimated x± from the 
measured Neel temperature. Instead, we could have de- 
termined X-L from the scaling measurements. In this way, 
all the parameters (x±, K y , and K z ) would have been ex- 
tracted from MQC measurements performed in the same 
regime of temperature, concentration, and field. How- 
ever, this would have required that we know beforehand 
the scaling of the excess spin (either from more sophis- 
ticated models, or from independent measurements). In 
the absence of such data, the use of the Neel temperature 
seems more reliable. 

As we have indicated earlier, A is of the order of 0.65 for 
fully loaded grains, which are therefore only close to the 
FM regime. However, A increases steadily with dimin- 
ishing loading, reaching a value of about 0.92 for a grain 
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with 1000 irons. Because of the magnitude of the ratio 
K z /K y , it turns out that, despite their quite large un- 
compensated moment, fully loaded grains are still within 
the range of validity of the small excess spin approach 
of Ref. O. However, this is no more true for partially 
loaded grains. 

One last remark about expression (|l6|) for the tunnel 
splitting. Strictly speaking, we see that it predicts a van- 
ishing splitting in the limit of small volumes. However, 
no hint of this appears in Fig. [|, which rather suggests 
that the splitting continues to increase. In fact, it is 
only below extremely small sizes (grain with a few tens 
of spins) that a sudden downturn occurs. Such sizes are 
not shown in Fig. |I| We note, however, that this is not a 
real effect, as the level splitting of small clusters has no 
reason to vanish, but rather signals a breakdown of the 
instanton approachEH. 



where J is the exchange coupling, k y and k z are the local 
anisotropics (felt by all spins). We are not aware of any 
molecular magnet for which all the magnetic anisotropics 
have been measured, hence we cannot give predictions 
for the tunnel splitting. However, many magnets have 
an excess spin of about the same order of magnitude as 
the sublattice spin, and since we expect k y and k z to be 
much smaller than J, this implies that KS/Safm 3> 1 
and that Awl. Therefore, they should be fully in the 
ferromagnetic regime. With respect to the results that 
would be obtained by neglecting the excess spin, notice 
that only the tunnel splitting is drastically different; the 
crossover temperature to the quantum regime will remain 
of the same order of magnitude. 

IV. DISCUSSION 



B. Molecular Magnets 

As we mentioned in the Introduction, racist molecu- 
lar magnets are ferrimagnets. For example^, the mag- 
netic properties of Mn 12 0i 2 (acetate) 16 point to a struc- 
ture where the spins of 8 Mn 3+ (5 = 2) are parallel to 
each other and antiparallel to those of 4 Mn 4+ (5 = 3/2). 
One can model this system as an AFM having two sub- 
lattices carrying each a spin 5 = 6, with an uncompen- 
sated spin 5 = 10 coupled to its Neel vector. LikewiseliS, 
[Mn(hfac) 2 NITPh] 6 , with 6 Mn 2+ (5 = 5/2) parallel to 
each other and antiparallel to 6 NITPh (5 — 1/2), can 
be viewed as an AFM with_s.ublattice spin 5 = 3 and 
excess spin 5 = 12. FmaUyB, [Vi 5 As 6 4 2(H 2 0)]- has 
15 spins 1/2 antiferromagnetically coupled, which in the 
AFM model gives a sublattice spin of 5 = 7/2 and an 
excess spin 5 = 1/2. 

To discuss the influence of the excess spin on the tunnel 
splitting in molecular magnets, it will prove convenient 
to express our results in terms of microscopic quantities. 
For this, we consider a perfect AFM containing N spins s 
(with coordination number z), whose Neel vector is cou- 
pled to an uncompensated moment formed by N e spins 
s. Our results read (again for k y <C k z ) 
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We have studied macroscopic quantum coherence in 
a simple model for a ferrimagnet, taken as an antifer- 
romagnet whose Neel vector is coupled to a large spin 
mimicking the uncompensated moment. 

We have obtained approximate instanton solutions for 
?i5/5afm large and K z 3> K y , where 5 is the excess 
spin, 5afm is the action of the instanton solutions in the 
absence of the excess spin, and K z > K y > are the mag- 
netic anisotropics. The expressions derived in this regime 
turn out to be also valid for ?i5/5afm small. Hence, our 
results cover all cases from FM to AFM. Known results 
are correctly recovered in both limits. Numerical simu- 
lations fully confirm the validity of our approximations, 
and this even for ratios K z /K y as small as 2. Using these 
solutions, we have calculated the MQC tunnel splitting 
by instanton methods, including the contribution from 
the fluctuation-determinant. 

Our results exhibit two surprising features. Intuitively, 
provided the uncompensated moment is much smaller 
than the sublattice magnetization (the tunneling entity) , 
one would expect only small modifications of the MQC 
resonance frequency. This criterion can be written as 
5 <C Ns, where 5 is the excess spin, N is the number of 
sites, and s is the local spin. In fact, we have shown that 
the correct criterion is much more drastic: The AFM 
regime requires that 5 <C Nsy/K z /J, where J is the 
exchange coupling. Since typically J is of the order of 
10 9 erg/cm 3 and K z of the order of 10 5 erg/cm 3 , there is 
a large reduction factor over the intuitive criterion. Even 
more surprisingly, for 5 3> Nsy/K z / J, one enters a FM 
regime, where the tunnel splitting is given by the expres- 
sions for a FM, but with the excess spin taking the place 
of the total spin. 

The quantity N s\J K z / J determines the critical mag- 
nitude of the excess spin beyond which one reaches 
the FM regime. We note that it can be rewritten as 
5afm y/K z /K y /%. For MQC to be observable, Safm/^ 
cannot be much larger than 10. Hence, the critical spin 
is quite small. This puts strong limitations on the pos- 
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sibility of observing pure AFM tunneling in mesoscopic 
grains. Ultimately, this shows that MQC is a very sen- 
sitive phenomenon and that a few additional spins can 
modify the tunneling frequency markedly. 

We have considered two applications. First, we 
have re-examined the experimental evidence for MQC 
in ferritin. Even though the magnetic moment of this 
molecule is small, it modifies the MQC resonance fre- 
quency already in the fully loaded grains used in the first 
experiments^, and has even stronger effects in the par- 
tially loaded grains studied laterE3. Our theory has sev- 
eral advantages over the original interpretation of the ex- 
periments. It does not make any hypothesis on the clas- 
sical attempt frequency ujq; instead, it predicts a value 
in fair agreement with recent measurements. Moreover, 
it reproduces the scaling of the resonance frequency as a 
function of the loading, whereas the original theory would 
not unless one was to arbitrarily set loq to a value much 
below the measured one. 

The value for the magnetic anisotropy extracted from 
the original MQC experiments was found to be— much 
smallerpthan the one derived from M6ssbauerE3 and 
SQUIDE3 measurements of the blocking temperature in 
the superparamagnetic regime. This was tentatively 
imputedt3 to differences in the fields and temperatures at 
which the various quantities are measured. With respect 
to the original interpretation, our theory gives a slightly 
larger value for K y , yet still much smaller than the one 
expected from the blocking temperature. 

At that stage, further testing of our theory requires ad- 
ditional data, either an independent measurement of the 
anisotropies in the MQC regime, or a measurement of the 
magnetic moment of the ferritin particles as a function 
of loading. 

The second application we have considered is to molec- 
ular magnets. Systems like these are usually neither pure 
FM nor AFM. Moreover, their small size renders them 
even more sensitive to an excess spin. Thus, our re- 
sults should be useful for a quantitative understanding 
of MQC in such systems. 
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FIG. 1. Instanton frequency and action as functions of 
TiS/Safm for K z /Ky — 30: numerical (diamonds) and an- 
alytical (line) results. 

FIG. 2. Instanton frequency and action as functions of 
HS/Safm for K z /K y — 5: numerical (diamonds) and ana- 
lytical (line) results. 

FIG. 3. Instanton frequency and action as functions of 
K z /K y for hS/S afm = 10: numerical (diamonds) and an- 
alytical (line) results. 

FIG. 4. MQC frequency as a function of the grain loading: 
experimental measurements (diamonds), and theoretical pre- 
dictions obtained with x± = 5 • 10 -5 emu/G cm 3 for a scaling 
exponent v = 1/3 (continuous line), v = 1/2 (dashed line), 
and with x± = 10~ 5 emu/G cm 3 , v = 1/3 (long dashed line). 
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